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We study linearized gravity in a six-dimensional Einstein-Maxwell model of warped braneworlds, 
where the extra dimensions are compactified by a magnetic flux. It is difficult to construct a 
strict codimension two braneworld with matter sources other than pure tension. To overcome this 
problem we replace the codimension two defect by an extended brane, with one spatial dimension 
compactified on a Kaluza-Klein circle. Our background is composed of a warped, axisymmetric 
bulk and one or two branes. We find that weak gravity sourced by arbitrary matter on the brane(s) 
is described by a four-dimensional scalar-tensor theory. We show, however, that the scalar mode 
£ — ' is suppressed at long distances and hence four-dimensional Einstein gravity is reproduced on the 

' brane. 
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I. INTRODUCTION 



Higher dimensional theories have been attracting much attention for many years. In the traditional Kaluza-Klein 
theories the compactification scale must be microscopic to guarantee effectively four-dimensional (4D) spacetime. In 
contrast to the Kaluza-Klein picture, recently proposed braneworld scenarios [l| allow for the presence of large (or 
even infinite) extra dimensions owing to localization of matter fields on the brane. Among various braneworld models, 
■ those with two extra dimensions have received particular interests from a phenomenological point of view. A naive 
\ reduction of six-dimensional (6D) gravity leads to Mpj = M 4 V @, where Mpi is the observed 4D Planck mass and 
M is the fundamental scale of gravity, and V is the volume of extra dimensions. For M ~ 1 TeV we roughly have 
' V 1 / 2 ~ 1 mm, which eliminates the hierarchy problem in the usual sense and at the same time offers us the possibility 
. to detect extra dimensions by table-top experiments on gravity. 

In the present paper we focus on the 6D brane model in which the internal two-dimensional space is compactified 
by a magnetic flux. (See, e.g., Ref. [3| and references therein for recent progresses in 6D braneworld models.) The 
basic motivation for considering such models is that extra dimensions, branes, and fluxes are important ingredients in 
string theory. Braneworld settings with flux compactification provide us natural toy models unifying aspects of string 
1 theory and cosmology, though such scenarios have not been explored much. Moreover, it has been suggested that the 
football (or rugby-ball) shaped codimension two braneworlds may help to resolve the cosmological constant problem 
•'■h , by a self-tuning mechanism (The mechanism has been criticized for several reasons Q.) 

The behavior of gravity in 6D braneworlds should be investigated carefully. In order for the brane scenario to 
be "realistic", it is necessary to reproduce 4D Einstein gravity on the brane at least on scales much larger than 
the compactification radius. Properties of gravity in a codimension two braneworld with flux compactification has 
been discussed so far in Ref. [f|. However, a codimension two brane suffers from the problem of the localization of 
matter. Namely, it is difficult to put energy-momentum tensor different from pure tension on the brane in Einstein 
gravity @, Q • Actually it has been demonstrated that massless relativistic particles can also be accommodated on a 
conical brane; the true problem here lies in the fact that the gravitational radius of a mass m is enhanced on the brane 
and hence gravity becomes nonlinear at distances much greater than a naive estimate [9( 1 . In any case the description 
of gravity on a conical defect is involved, which significantly restricts the analysis of gravitational aspects of 6D 
braneworlds. One way to evade this generic problem is introducing codimension one brane with one extra dimension 
compactified on a circle (which we call an extended brane) instead of the original codimension two defect PH. (See 
Refs. d, [ll|, EH EH E3, EE Hl| for other ways out.) Since the model has a Kaluza-Klein spatial direction in addition 
to one large extra dimension, this construction is thought of as a "hybrid" braneworld scenario p7| . (Very recently 
this approach has been employed to regularize a conical singularity in a different context [l8j.) 
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Along this line, the authors of Ref. [10( were the first to study the behavior of gravity in the presence of matter 
sources on the brane. In their model the internal space is football shaped (i.e., not warped) and codimension two 
branes at the poles are replaced by the extended branes wrapped around the axis of the football. The perturbation 
analysis shows that standard 4D gravity is reproduced on the extended brane [l(| • The system is stable under the 
most generic axisymmetric perturbations [19]. In Rcf. [2Cj the same regularization technique is used for more general 
brane models with warped internal space both in Einstein-Maxwell theory Q and in 6D supergravity [U [22j • In this 
paper we extend the work of [lOj and investigate linearized gravity sourced by arbitrary matter on extended branes in 
the Einstein-Maxwell model of the warped braneworld [2(| ■ As in [l(| , we basically follow the well-developed analysis 
of perturbations [23| in the Randall-Sundrum braneworld (24j. This model will serve as a simple playground to study 
gravitational aspects of braneworld scenarios in the context of string flux compactifications. 

The paper is organized as follows. In Sec II wc describe our background model. Then wc derive perturbation 
equations and boundary conditions in Sec. III. In Sec. IV we discuss the behavior of linearized gravity on the brane. 
Finally we draw our conclusions in Sec. V. 



II. 6D WARPED FLUX COMPACTIFICATION AND EXTENDED BRANES 



A. Bulk solution with extended branes 



The 6D background model we will discuss was originally used to compactify an internal two-dimensional space in 
Einstein-Maxwell theory [25|. We consider the system in the braneworld context, and basically follow the construction 
of [10, 20]. The bulk has a structure of Minkowski spacetime times a two-dimensional compact space, and consists of 
M + 1 regions (Mi) separated by M extended branes (Si)- O ur action is 



(2.1) 



where the 7-th bulk and i-th brane actions are given, respectively, by 



Si = / d b x^g 
J Mi 

Si = — d 5 xy/~^q 

J-Ei 



A/ 4 



MN 



(2.2) 
(2.3) 



Here the capital Latin indices {M, N, ■ ■ ■ } and the Greek indices with hat {//, £',■■■} are used for tensors defined in 
the 6D bulk and on the 5D brane(s), respectively. We have the U(l) gauge field Am with field strength Tmn = 
9m An — 9m Am in the bulk, which is coupled to the brane scalar fields Cj. On each brane Oi arises as a Goldstone 
mode of a Higgs field, and the absolute value of the vacuum expectation value of the Higgs field is given by Vi/y/2. We 
also include tension Ai in the brane actions 2 . In this model the cosmological constant 1/(2X|) is different at different 
regions in the bulk. 

The background solution is described by the metric 



gMNdx dx = a (w)rj t j lV dx^dx 1 ' + L I 



dw 

7R 



Pif(wW 



where 



and the field strength 



1 



5(1 -a) 2 



-a 2 (w) + 



a 3 (l-a 5 ) 1 



1 — a 3 a 3 (w) 



a 6 (w) 



(2.4) 



(2.5) 



fiM 2 (3iL! 



! 3a 3 (l-a 5 ) 
5(1 -a 3 ) ' 



(2.6) 



2 The brane action of the form ( 12.:{[ ) may not be unique but rather is the simplest one. See Ref. [2ll for brane actions. 
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The warp factor is given by 



a(w) = - [(1 - a)w + 1 + a] , 



(2.7) 



and a (0 < a < 1) characterizes warping of the bulk. The function f(w) has two positive roots w = 1 (a = 1) and to = 
— 1 (a = a), and hence we will consider the space — l<u><l(a<a<l). The original construction of the braneworld 
model of this type is given by Mukohyama et al. in Ref. [6J, in which codimension two defects are considered. The 
stability analysis against classical perturbations has been performed in their subsequent papers [2(||27j]. In Appendix El 
we replicate the detailed derivation of the bulk solution [6l. [20l|. 

We put brane(s) at w — Wi. Then the continuity of the {4><j>) component of the induced metric implies that PiLj is 
continuous across the brane. Therefore we write 



PiLi 



(2.8) 



The points w — ±1 can be regarded as poles. We assume that </> has period 2tt. In order to avoid conical singularities 
at the two poles, we impose 



20(l-a)(l-a 3 ) 
5 - 8a 3 + 3a 8 : 



P- = 



20(1 -aT 1 )^- a" 3 ) 
5 - 8a" 3 + 3a" 8 : 



(2.9) 



where (3± denotes the value of /3j in the region that contains the pole w = ±1. It is easy to see that f3 + > 1 > /3_ > 0. 
Now we check that the above bulk solution coincides with the model of [X 01 ] in the unwarped limit a — > 1. Indeed, 
we have / — > 1 — w 2 , T w ^ — > M 2 PjLj, and (3± — > 1 as a — > 1. With the coordinate transformation w — sin 9 we can 
reproduce the bulk solution of 10]. In this unwarped limit the background has a symmetry across the equator 
(w = 0) before introducing branes. 

Let us look at the quantization condition for the flux. Since T wi p — A'^, where the prime stands for a derivative 
with respect to w, we have 



2fi£M 2 
3(1 -a) 
2(j,£M 2 
3(1 - a) 



1 

a 6 

1 ] 



for w > w e — e, 

for W < W e + £, 



(2.10) 
(2.11) 



where w e is arbitrary and e > 0. The integration constants here are chosen so that = at the poles. We see that 

in the overlapping region w e — e < w < w e + s, A</, is doubly defined. This means that A^ and A^f are related via 
the following gauge transformation: 



A 



(S) 



A N) + (^M 2 <t>) - 



'(l-a 3 )(l-a 5 ) 
15a 3 (l-a) 2 



(2.12) 



Noting that Oi is the phase of the brane Higgs field with charge e under the t/(l) symmetry, we arrive at the following 
quantization condition: 



2eiiJM 2 = N, N = 0,1,2,- •• 



(2.13) 



The junction equations at a brane are given as follows. They are derived from (i) the continuity of the induced 
metric, (ii) the Israel conditions, and (iii) the junction conditions for the Maxwell field. The first one was already 
used to relate /?/ with Lj in Eq. (|2.8[) . The Israel conditions are written as 



' M 4 ' 



where S^o is the energy-momentum tensor on the brane and [[-P 1 ]]^ := lim £ ^o (^l^i+e — ^| tB ; 
to compute the extrinsic curvature. The energy-momentum tensor on the brane is given by 



(2.14) 

s ). It is straightforward 



c(«) _ 



& A, A, — 



Aj + -^qf* {d^ - eA 4 



A,: - -j-q^ (8^ - eA+Y 



(2.15) 
(2.16) 
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FIG. 1: Single brane and two-brane models. 



Then the Israel conditions are rearranged to give 



f 1/2 (f 



l \2f a 
4/ 2 a 



The junction conditions for the Maxwell field can be written as 

[[n N ^nm]]^ i = - ev 2 (d M <Ji - eAu)^. , 
where n N is the unit normal to the brane. For the background, only the <j> component is relevant: 



From Eqs. (f2TT7|) and (|2~20| we obtain 



8 fl/2 [ / a 

2/ a 



(2.17) 
(2.18) 

(2.19) 
(2.20) 

(2.21) 



Since /'/(2/) — a' /a = has a root w c (—1 < w c < 0), we see that [[/Sr]]^ > for Wi > w c and [[/Jr]]-^ < for 
Wi < w c . 

The equation of motion for Oi implies (d^di ~ eA^) = 0, leading to <Ji — rii(f>, where rii is an integer. Using 
Eqs. (f2~T7]) and (|2~20]) we obtain [H 



N 



I -a 3 



5(1 -a 8 ) 3 
8(1 - ofi) 1 



(2.22) 



where c, = 1 for Wi > w e and Ci = a 3 for Wi < w e [see Eqs. (|2.10p and (|2.11[) ]. In the limit a — > 1, we can see that 
rii = —N/2 for Wi > w e and rii = +N/2 for v)i < w e [Tfjj . 
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FIG. 2: The position of the brane in the single brane model, 
brane below this line. 



The dashed line shows the value w c and hence one cannot put a 



B. Single brane model 

Let us consider a model with a single brane at w — w. The two regions separated by the brane have the curvature 
scales L + = £/ ' j3 + and L_ = £/(3-, respectively. The configuration of the model is described in the left side of Fig. [T] 
For given parameters a and £(ev) 2 , Eq. (|2.2ip determines the position of the brane w, as is shown in Fig. Note 
that one cannot put the brane in the region w < w c because of the condition /3 + > /3_ . 



Two-brane model 



We can also consider a two-brane model, which is nothing but the situation considered in Ref. [20(. In this case, 
the bulk is divided into three parts by two branes at w — w+ and w — w_ . A schematic picture is shown in the 
right side of Fig. Q] The three bulk regions are characterized by the curvature scales L+ = £//3+, Lq := £//3q, and 
L_ = £/ '/?_, respectively, from north to south. If w c > w+ > W-, it follows from Eq. (12.21j) that f}+ < (3q < f3~. This 
is incompatible with the fact /3+ > /3_, and so we impose vi+ > w c . Consequently, we have (3+ > (3q. However, both 
w- > w c (0q > j9— ) and ui- < w c (0n < /3_) are possible. 

The background configuration of [1 01 ] can be recovered by taking the unwarped limit (cv — ► 1) and assuming at 
the same time that w+ = u_. As is mentioned in the above, the unwarped bulk automatically has the .^-symmetry 
across the equator at w = 0. The second condition requires that the branes are located at Z2-symmetric positions, 
i.e., W- = —w + . Furthermore, in Ref. [Io| the ^-symmetry is imposed also {or linear perturbations. We stress that 
the general warped background (a / 1) has no such symmetry. Our perturbation analysis in the rest of the paper 
does include the unwarped model as the limiting case, but we will not restrict ourselves to the ^-symmetric brane 
configuration nor perturbations in the limit a — > 1. We will comment on this point in more detail in Sec. HVBl 



III. LINEAR PERTURBATIONS 



We now consider linear perturbations from the background solution described in the previous section. We write 
the perturbed metric in an arbitrary gauge as 



(9 



MN 



o~9mn) dx dx = a 2 [{l + 2^)ri^ + 2E^ + E^ + E l ,^ + h^]dx^dx v 
+2 {B Wtll + B WI1 ) dwdx» + 2 (fl „ + B M ) d^dx" 



(1 + 2W) 



dw 2 



2C'Pjfdwd(l) + (1 + 23>)/3 2 I fd(f) 2 



(3.1) 



where the perturbations are split into scalar, vector, and tensor modes under the Lorentz group in the external 
spacetime. The vector modes E^, B wfi , and B^ are transverse: E = 0, ■ • ■ . The tensor mode is transverse and 
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traceless: h' v = = h A Similarly, we can write the perturbations of the gauge field 5 Am as 

SAf, = df+A + A^, SA W = A w , SA^ = A$, (3.2) 

where A^ is a vector mode and the other three are scalar modes. The scalar part of the perturbed field strength 
8T MN is given by 

wrj> = A^, SJ- ncf, = A^^, 8J- n W = F w ^, (3-3) 

where F w :— A w — A' . The vector part of the perturbed field strength is 

^ M = i;, 5F tlu = d fl A v -d u A IM . (3.4) 

We shall discuss tensor and scalar perturbations in the following subsections. We show in Appendix IBI that there are 
no vector-type perturbations of importance. 

A. Tensor mode 

Tensor perturbations are invariant under an infinitesimal coordinate transformation. The 6D equation of motion 
for tensor perturbations is 

(a 4 A^)' + a 2 £?nV=0, (3-5) 
where □ := rf iV d ii d v . Setting □ = we obtain the zero-mode solution 

/W (jfljj 
irr (3-6) 

where k 1 and B 1 are integration constants. Requiring the regularity of h^ v and h' v at the poles [H, |27|, we see 
that B 1 ^ vanishes in the region that contains the pole. Thus, in the single brane model the zero mode is given by 
/i M „ = k^ix) everywhere. The source-free Israel condition reads [[/Jj/i^,,]]™; = 0. This means that, also in the two- 
brane model, the zero mode is written as /i M „ = A^(x) everywhere. In the next section, we will further discuss the 
tensor mode with the introduction of matter sources on the brane (s). 



B. Scalar modes 



This section focuses on the master equations and boundary conditions governing the dynamics of scalar-type 
variables. The details of the computations of the Einstein and Maxwell equations are given in Appendix [Cl 

As shown in the case of the codimension two model [H,[13j the perturbation equations in the longitudinal gauge 3 , 
(|C20j) - (|C24|) and (|C27j) . can be solved using two master variables defined by 

fil := -$ - 3#, Q 2 ■= (3.7) 
With a straightforward computation Eqs. (|C20|l - (|C24"|> and (fC27)) can be recast into a set of coupled equations 



L 
f 



5— mi - - (n 



L 



a 2 f 



4-n' 



L 



f-D(ni+20 2 



0, 
0. 



a ~ 2a 2 /" 

In terms of these master variables, the remaining metric and gauge field perturbations are given by 

W = Qi+fi 2 , 

H 



£M 2 a A 



^ = f 



n[ + L (fix + 2n 2 ) + 2-n! 

/ « 



(3.8) 
(3.9) 

(3.10) 
(3.11) 



3 Scalar-type variables with tilde denote perturbations in the longitudinal gauge, which is explained in Appendix [Cl 
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We have three more variables, C, A, and A w . As explained in Appendix [Dl we do not have to take care of these modes 
in any case. 

We are interested in the zero-mode sector because of its importance in describing the long range gravitational force. 
Setting □ = in Eqs. (|3.8p and (|3.9[) . we obtain analytic solutions for the zero modes: 



1 



5(1 -a) 2 / 



pi(x) ( qj(x) , r ^ 2 Avi{x) 



ui(x)a 



ui(x) 



vi{x) 



(3.12) 
(3.13) 



where the integration constants pi,qi,uj, vi are to be determined by the boundary conditions. In the unwarped a = 1 
case, general solutions for the zero modes are 4 



1 



Pi(x) + qi(x)w + U[{x)w 2 + —Vi(x)w 



1-w 2 

^2 = Uj(x) + Vj(x)lV. 

Since / = at the poles and fii should be regular there, we impose 
We also require that A$ = at the poles, i.e., 

[(/fii)' + 2/ / n a ] U)=±1 = o. 



(3.14) 
(3.15) 

(3.16) 
(3.17) 



A more rigorous argument on the boundary conditions at the poles is found in Ref. [27j. 

We now turn to the boundary conditions at the branes. The continuity of the induced metric and the gauge field 
imposes 



a 

V ' 
$ + — C 

2/. 



= o, 
= o, 
- o, 



(3.18) 

(3.19) 
(3.20) 



where we simply write the bending mode of each brane as = £. 

The (fiu) component of the Israel conditions, including the transverse and traceless perturbation, is 

Li „ n/ , / 1/2 



f 1 ' 



J- LIV 



(3.21) 



where T$ is the matter energy- momentum tensor on the brane labeled by i. The four-dimensional trace of Eq. (|3.21|) 



reduces to 



Lr 1 



The 



/V2 a 2 

component of the Israel conditions gives 

7 1/2 



IT. 



Mi) 



3 M 4 



Li 



_ 4 *' + A-W + ( L - -) $ + 4W + Y } 
a \2f a J a z f J 



T, 



M 4 



(3.22) 



(3.23) 



4 Eqs. I|3.12jl and l|3.13| l do not coincide with Eqs. H3.14H and (|3.15jl in the limit of a — > 1. However, by replacing the coefficients ui and vj 
as uj — > uj + [(5 — 3a)/3(l — a)]' u r and — > —[2/3(1 — ct)]v[ we obtain an expression for Q2 that reduces to Eq. 113.151 1 in the unwarped 
limit. We can obtain such an expression for f2i by replacing pi and qi in Eq. l|3.12j l with appropriate combinations of pj , qi ,uj , and 
vj. The resultant expression is unnecessarily complicated. 
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where we defined 



1 "-is- 7) (f f -4)< + ^f< + 7m^^ 



The component of the Maxwell junction conditions gives 



= 0. 



(3.24) 



(3.25) 



IV. LINEARIZED GRAVITY ON THE BRANE 



A. The zero-mode truncation 



Following the approach of [23[ we shall see the behavior of weak gravity created by matter sources on the brane(s). 
Rearranging the Israel condition (|3 . 2 1[) we obtain 



where we collected together the matter sources and the brane bending scalars in the right hand side: 



cS (i) — 2a 2 \ l J±l. 



T- 



A(i) 



(4.1) 



(4.2) 



with aj := a(wt) and ft := f(wi). We can now put the equation of motion (|3 . 5[) and the boundary condition (|4.1 
into a single equation with source terms as 



i 

where we defined an operator 

£>V := [a%]' + cpe^Uhp,. 
Then using the (retarded) Green function defined by 

OGr(x, w; x', w') (a; — x')S(w — w'), 

we can solve Eq. (|4.3[) : 

h fJlV (x,w) = - J d 4 x' G R (x,w;x',Wi)S$. 

i 

The Green function is explicitly given by 



Gr(x, w; x , w 



dik ik-(x-x') ST^ lpn{w)j> n (w') 

(2tt) 4 ^ ml + k 2 - (to + ie) 2 ' 

where tp n (w) are a complete set of eigenfunctions of 

and are normalized as 

/ ipn(w)ip n '(w)a 2 (w)PY 1 dw — S nn > ■ 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



Here we are interested, for example, in gravity acting on isolated sources separated by a distance much larger than 
the bulk length scale I. In such situations, we may expect that the effect of Kaluza-Klein modes is suppressed so that 
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we are justified to truncate the Green function by keeping only the zero-mode contribution. As already shown in the 
previous section the zero- mode solution is constant, ipo — £~ 2 , where 



I3j i a 1 {w)dw 



1/2 



Note that we basically have ~ 0{t). The zero-mode truncation of the Green function yields 

h // m ) 4- 

where we have separated /i M „ into the "matter" and "brane bending" contributions as 



2f 



T 



Hi) 



£ 2 M 4 
2£ 2 



2 ^iiu i ' 



G 



(4.10) 
(4.11) 

(4.12) 
(4.13) 



In order to discuss the behavior of gravity on the brane, we compute the linearized 4D Ricci tensor for the induced 
metric 



q® = of [(1 + 2*) r illv + V]L 4 . 
where * is a metric perturbation in the Gaussian-normal gauge. We can write the Ricci tensor as 



where 
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(i) 



2* + -n/i (c) 



(4.14) 



(4.15) 



(4.16) 



Using the trace of the Israel condition (|3.22p . we find that the first two terms give 



T 



Hi) 



T 



Hi) 



2 ^jui> 

where T Ml/ is the energy-momentum tensor integrated along the direction, 

Jo 

and the 4D gravitational coupling at each brane is defined as 



(4.17) 



87rG (l) 



a 2 {wi) 
2it£ 2 M i ' 



(4.18) 



(4.19) 



Now it is clear how the 4D tensor structure can be recovered by the first two terms in Eq. (|4. 1 5[) : the brane bending 
plays a crucial role. The same mechanism works in the Randall- Sundrum model [23| and in the unwarped, Z2- 
symmetric model of [h| . What would make brane gravity different from 4D Einstein is the ip^ term. Apparently, 
gravity on the brane is described by a scalar-tensor theory due to this scalar mode. However, as we shall show in the 
next subsection, the effect of <p^ term can be ignored as far as long range gravity is concerned. 
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FIG. 3: Plot for A(w, a). The location of the brane is restricted in the region w > w c 



B. The <p (i) mode 



We will now evaluate the impact of <p(*' on brane gravity and make clear the physical meaning of this mode. We 
start with combining Eqs. (|3.22|) and (|3.23|) to obtain 



fa I -W + 4-W + ( - - ) $ + Y 
a \2j a 



(4.20) 



where 



(i) ._ if Hi) _ rp Hi) 



(4.21) 



and T^ 1 ^ is defined similarly to T^J . This combination of the energy-momentum tensor is the key quantity. To catch 

the meaning of (p^ l \ the following two quantities are important: the perturbed volume of the internal space SV and 
the perturbed circumference of each brane SC^ 1 '. We can compute the volume of the internal space by performing the 
integral 



~Ai) 



Vq + SV = / g ww g^dwd<j), 



(4.22) 



yielding V = 2tt£ 2 J j3j dw and 



SV = -2tt£ 2 + 2 J pj^dv^j 



The circumference of the brane is given by 



and so we have Cq^ = 2tt£ ft'* an d 



fl/2 



4° + SC® = 



6C® = 2ir£fl /2 $1 



(4.23) 



(4.24) 



(4.25) 



Let us first consider a single brane model. In this case there are 4x2 = 8 scalar integration constants and 2 brane 
bending modes. Using 4 boundary conditions at the brane [Eqs. (|3.18|) - (|3.20|) and (|3.25p ] and 2x2 = 4 regularity 
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FIG. 4: A (+) as a function of w+ and a. The other parameters are fixed as W- 
restricted from the condition w+ > m&x{w-,w c }. 



= —0.8 and /3o = 0.8. The domain of w+ is 



conditions at the two poles [Eqs. (|3.16[) and (|3.17[) ]. we can express 8 of 10 variables in terms of 2 variables, say u_(x) 
and u_(x). Then, substituting the result into Eqs. (|4. 16|) and ()4.20p we can express ip and f in terms of U-(x) and 
V- (x) . This procedure reveals the relation 



a 2 (w) 



8ttG? A(w,a), 



(4.26) 



where A(w, a) is a regular function and it is dimcnsionless. Since the explicit form of A(ui, a) is quite messy, we show 
a plot in Fig. [3] instead of writing down a lengthy equation. Now we have 



Ruv — 8ttG 



T, 



-^q^ )-A(w,a) 



a 2 (w) 



1 



(4.27) 



From this it is easy to see that the second term is negligible on scales much larger than (~ t) . 

The scalar mode tp decouples from gravity on the brane as the brane shrinks to the pole (w — > +1), since A — ► 
in this limit. This is in agreement with the result of (lp| . Note, however, that in fact one cannot bring the brane to 
the pole for a finite value of v because of Eq. (|2.21[) in the single brane model. 

We can also express SV and SC in terms of u_(x) and u_(x). With this we find 



tp oc 5V oc SC. 



(4.28) 



This means that the perturbation of the volume of the internal space and that of the circumference of the brane are 
not independent, and tp corresponds to this mode. 

The situation for the two-brane model is similar to that of the above single brane case. There are 4 x 3 = 12 
integration constants and 2x2 = 4 brane bending scalars, while we have 4x2 = 8 boundary conditions at the branes 
and 4 regularity conditions at the poles. Thus we can express 12 of 16 variables in terms of 4 variables, say uo(x), 
i>o(x), u_(x), and u_(x). Then Eqs. (|4. 1 6[) and (|4.20|) allow us to express p^ + \ f( + ', and f in terms of these four 
variables. Interestingly, with this procedure we can show that all the integration constants in p^ are encoded solely 
into f^' as 



(+) = 1±. 8jtG (+) 
a± 



F(+) A (+) +f (-) A (-) 



(4.29) 



where A^' are regular functions of w+, W-, a, and (3q. Thus we arrive at the same conclusion as the above: the effect 
of the mode tp^ can be safely neglected on scales much larger than on the brane. Again, the explicit form of 
A^) is quite involved and so we just illustrate their behavior in Figs. H]-[7] It turns out that A( + ) — ► as the brane 
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FIG. 5: A' ' as a function of w+ and a. The other parameters are fixed as W- = —0.8 and /3o = 0.8. The domain of w+ is 
restricted from the condition w+ > max{u)_,w c }. 




FIG. 6: A' + ' as a function of W- and a. The other parameters are fixed as w+ — 0.8 and f3o = 0.6. 



at w — w+ shrinks to the north pole, w+ — > +1. Similarly, we have A^ - -* — > as the brane at w = w_ shrinks to the 
south pole, W- — ► — 1. Thus, the <^( + ) mode decouples from gravity on the brane when one brings both branes to the 
poles. 

It is straightforward to express the perturbations of the circumferences of the two branes, 5C^ and 6C(~', and 
the volume fluctuation 5V, in terms of uq(x), vq(x), U-(x), and V-(x). Comparing the result with f^' expressed in 
terms of these four variables, we find that 6C& and 5V can in fact be written as linear combinations of r' •': 

ae(±) = 4 ±) f(+) + c ( ±) f(-', (4.30) 

SV = dif< + > +d 2 ? i --\ (4.31) 

Namely, two of the three fluctuation modes 5C^ ' and SV are linearly independent. The mode can be interpreted 
as a combination of two of these fluctuation modes. 

Finally, let us take a brief look at the unwarped a = 1 model, focusing on the simplest background configuration 
with u + = V—, i.e., the ^-symmetric configuration of the branes {w— = — ?Z>+). In this case it is instructive to 
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FIG. 7: A' ' as a function of W- and a. The other parameters are fixed as w+ = 0.8 and f3a = 0.6. 



rearrange Eq. (|4.29p as 
The coefficients are given by 

A« - A(") 



A(+) + A(_) r fM + f (-)) + A(+) - A ^ ( f M _ f i-) 



(1 - w+) [9/3 + (10 - 7(3)w+ - (11 + [3)w+(l + «;+)] 
12 [/?(1 - tD+) + w+] 
w+(l - w+)(10 + w+ + g+) 
2 [10/3 + 9(1 - /?)«>+ + (1 - [3)w%] ' 



(4.32) 

(4.33) 
(4.34) 



The restriction to .^-symmetric perturbations pro jects out the second term, and then the result is in agreement 
with [10| . In other words, Eq. (|4.32|) generalizes [10| to non-Z 2 -symmetric perturbations. 



CONCLUSIONS 



In this paper we have studied linearized gravity on a brane in a 6D Einstein-Maxwell system. In the present 
braneworld model two extra dimensions are compactified by a magnetic flux, and hence the model is thought of as a 
toy model of string flux compactifications. Introducing codimension one branes (extended branes) with one spatial 
dimension compactified on a Kaluza-Klein circle instead of strict codimension two defects, one can put arbitrary 
energy-momentum tensor on branes [l(J [l!| [2(|. With this we can first discuss the behavior of gravity sourced 
by matter on the branes. We have considered an axisymmetric, warped background with one or two branes, and 
examined axisymmetric perturbations in this system, generalizing the previous analysis [10( to the warped model 
of • We have focused on the zero- mode sector of perturbations as the effect of Kaluza-Klein modes is expected to 
be subdominant at long distances. 

It is a well-known fact that the contribution from the brane bending mode is crucial for recovering the 4D tensor 
structure in the Randall-Sundrum braneworld [23j . We found that the same mechanism works in the present class 
of 6D warped braneworlds as well, which has been shown only for the unwarped Z2-symmetric model in [10| , There 
appears another scalar mode associated to the fluctuations of the circumferences of the branes, so that gravity is 
described by a scalar-tensor type theory. However, we have shown that, as in the unwarped case, the scalar mode can 
be neglected at distances much larger than the scale of flux compactification, ~ t. Therefore, 4D Einstein gravity 
is reproduced on an extended brane in 6D warped flux compactifications. 

On scales about I the effect of Kaluza-Klein modes will not be negligible in addition to the correction from the 
above scalar mode. The precise evaluation of the corrections from the Kaluza-Klein modes is important for testing 
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the model against gravitational experiments and observations, which is one of our remaining issues. Note, however, 
that the scale of the circumference of the brane is given by £f 1 ^ 2 (wi), which is much smaller than the scale of the 
internal space £ provided that the location of the brane is sufficiently close to a pole. Consequently, the effect of the 
modes which are inhomogeneous along the <fi direction will become significant on much smaller scales than £. Note 
also that even though we have one microscopic Kaluza-Klein direction, the 4D Planck scale is given by Eq. (|4.19p 
[not by Mpi <~ £ ■ £f 1 / 2 (wi)M 4 ]. Thus it is possible to have the scale £ of order submillimeter while addressing the 
hierarchy problem. 

Finally we would like to make several remarks on possible extensions of our current work. First, the stability analysis 
for the present model remains to be done. It would be also interesting to explore cosmology and nonlinear brane 
gravity on an extend brane, as we are now able to put arbitrary matter on the brane. The cosmic expansion can be 
described by a moving brane in a warped bulk, with the scale factor mimicked by the warp factor (28j . However, in the 
present model the warp factor is bounded above and hence it seems difficult to get an ever-expanding brane universe 
following this line. Note that a de Sitter universe on an extended brane has been constructed straightforwardly for 
the unwarped bulk jT^]. (See also Ref. @.) Going beyond the Einstein-Maxwell model, we also plan to study aspects 
of gravity in the supersymmetric version of the present brane model [2(| HH, [HJ . 
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APPENDIX A: DERIVATION OF THE BACKGROUND METRIC 



We consider the 6D Einstein-Maxwell system described by the action 

' M 4 ( _ 1 \ 1 



S = / d b Xy/=g~ 



The field equations derived from this action are 



2 \ n -v)-\ T2 



(Al) 



TImn - -9mnT^ - ~^jp9MN + -^j (^Fm L F n l - -9mnF 2 ) ■ t A 2 ) 

We consider a solution obtained by a double Wick rotation from the Reissner-Nordstrom solution [6j, |26j, |22j . The 
metric can be written as 

dn 2 

g MN dx M dx N = p 2 r,^dx»dx v + J- + f{p)d4> 2 , 

fyp) 

f! \ - D P 2 Q2 1 (A^ 

nP) ~ P 3 20L 2 12MV 6 ' 1 j 



and the field strength is given by 



We assume that the metric function f(p) has two positive roots p+ and p- and f(p) > for p_ < p < p + . The 
integration constants D and Q are conveniently parameterized by a := p-/p+(< 1) as [20| 

^l_a« n2 _ p\M 4 3a 3 (l-a 5 ) 

20L 2 l- a 3' L 2 5(1 -a 3 ) ' [ ' 

Introducing a new coordinate w defined by 

a(w):=l[(l-a)w+l + a} = -?-, (A6) 
2 p+ 
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we can rewrite the metric and the field strength as 

gMN<lx M dx N — a 2 p 2 + rj liv dx^dx 1 ' + 



2L 



(1-a) 



1 M 2 / 3a 3 (l-a 5 ) 
* ~ ^TV 5(1 -a 3 ) ' 



where 



1 



I 2 , (A7) 
(A8) 

(A9) 
(A10) 

We would like to remark that in the above the Minkowski metric rj^u can be replaced by any Ricci flat metric g)w ■ 



1-a 8 1 a 3 (l-a 5 ) 1 
1 — a 3 a 3 



1-a 3 a b 



5(1 - a) 2 

Note here that / depends only on the parameter a. The rescaling of the coordinates 
leads to the background solution in the main text. 



APPENDIX B: VECTOR MODES 



Under a vector gauge transformation 
the metric perturbations transform as 





-> x" - 










B — 













From Eqs. (|B2|) wc find the two gauge invariant variables 



Vft :— B wlx ~ a 2 E^ 



and The perturbed gauge field A^ is also gauge invariant. 
The linearized Einstein equations give 



V' + l-Vp + ijVp = 0, 



0, 



a' , fa'\ 2 L 2 , I 2u j. 

From the \x component of the linearized Maxwell equations we obtain 

M 2 p 



The {pv) component of the Israel conditions gives 



Ml 

2£ 



h (V^ + V v ,, 



LpB„ = 0. 



M 4 ^ v ' 



where a, := a(wi) and fi := f(wi). From the (fJ,<p) component of the Israel conditions we have 



1 , a' \ /Va 
a 



ev 2 (d<f,a - eA<j>) A p 



(Bl) 
(B2) 
(B3) 

(B4) 
(B5) 
(B6) 

(B7) 

(B8) 
(B9) 
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The Maxwell junction condition is 



e(ev t ) 2 A, 



(BIO) 



Equations (|B6p and (|B7[) with (|B9|) and (|B10|) govern the two vector variables and However, they do not 
couple to matter on the brane via the junction conditions. For this reason we will not take care of these modes. 
Equation (|B5[) indicates that only the zero mode is present for V^. Equation (|B4p is then solved to give 



a 2 f ' 



(Bll) 



where is an integration constant. In order to ensure the regularity at the two poles, we impose = in the region 
that contains the pole. In the single brane model, this immediately leads to Vn = everywhere. In the absence of 

vector-type matter sources, the source-free Israel condition, Eq. (|B8[) with Tf$ — 0, forces to vanish everywhere 
in the two-brane model. 



APPENDIX C: PERTURBATION EQUATIONS AND GAUGES FOR SCALAR MODES 

1. Gauges 

Under a scalar gauge transformation 



the metric perturbations transform as 



- 




w - 


- w + C, 


cj> - 


- 4> + ^, 



(CI) 



The perturbed gauge field transforms as 



and so 



* -> 


* - 








a 


E -> 




£, 


B w — > 






B -> 


B - 




VF -> 


W - 


t 10 ' _J_ ' 



r 2 



f 

<f> y <f> — 



a - ^-r^i, 



F F _|_ /I' f4 



(C2) 



(C3) 



(C4) 



One of the gauge we employ in the present paper is an analogue to the longitudinal gauge, which is defined by 
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E = B w = B = 0. Hence we can define metric perturbations in the longitudinal gauge as 





:= V — 






vv 


. — vv — 




/; X 


$ 


■= $- 


f 

2Lj X > 




C 


:= C- 







where wc defined a convenient variable X 
gauge can be defined as 



(C5) 
(C6) 
(C7) 
(C8) 

B w — a 2 E' . Similarly, the gauge field perturbations in the longitudinal 



A := A - As 
Aa ■= A 



B 



A'^jjX. 



(C9) 
(CIO) 
(Cll) 



Scalar-type variables with tilde will denote perturbations in the longitudinal gauge. 

The 6D perturbation equations in the longitudinal gauge reduce to a system of two (coupled) master equations. 
For this reason the longitudinal gauge is useful for solving the bulk equations. In this gauge, however, the positions 
of the branes are perturbed in general: 



10; 



(0 



w 



>+C W (*). 



(C12) 



Strictly speaking, we have two different brane bending modes for each one brane, C™>m>; an d Cw<wi > as no -^-symmetry 
is assumed at the branes. 

The Gaussian-normal gauge, defined by 



B n 



W = C* = and w 



(i) 



(C13) 



is convenient when imposing the boundary conditions at the brane as the brane position is not perturbed. We 
can employ this gauge in the neighborhood of each brane. Hereafter, scalar-type variables with bar will denote 
perturbations in the Gaussian-normal gauge. Starting from the longitudinal gauge, the Gaussian-normal gauge is 
realized by a coordinate transformation x M = x M + £ M such that 



w 



c-e , 
C (i) + c 



(C14) 

(CIS) 

(C16) 
(C17) 



We can fix the residual gauge freedom by imposing ^| tBi = £^\wi — 0. With this, perturbative matter sources localized 
on the brane are invariant when going from the longitudinal gauge to the Gaussian-normal gauge. 

The above gauge transformation relates the longitudinal gauge perturbations with the perturbations of the induced 
metric and gauge field perturbations on the brane. For example, we find 



*k = 



$1 



t a a 

* + — c 

a 



(i) 



f 

2P 



(C18) 
(C19) 



We also have A\ 



Aim. and F. 
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2. 6D perturbation equations for scalar modes 

We work in the longitudinal gauge for the analysis of the bulk perturbations. The linearized Einstein equations 
give 



-f : : (<!• -2/ 

4/ 



4 + 6-)*' 

/ a 



a 8 £M 2 a 4 



v l' " - ( £ + 5-) *' - -W" + - ( 4 + 3- ) - 
2/ a / a a \ ,/ a ' 



L 2 ( - ~\ - f ~ 
a 2 V / a 8 



lM 2 a^ 



*" + < 27 + 4) + s (£ + 7) 41 + 3 (y + 4 1 5 " 



^ + 6-4 + 6^)(* -in 



^(n + $) 

a 8 



1 + = 



/i 



4' 



and 



2* + $ + n = o, 



□c = 0, 



f" + '>[£■ + c -. 



2/i 



£M 2 a 4 

From the linearized Maxwell equations 8mS [y/—gF MN ] — we get 

' a 2 L 2 



a A'. + fiM 2 £(W - $ - W) 



+ — p-n^ - 0, 



and 



0, 
0. 



(C20) 

(C21) 
(C22) 



(C23) 
(C24) 

(C25) 
(C26) 

(C27) 

(C28) 
(C29) 



On each brane we have the perturbed equation of motion for <7j field: \_\f- ~q_q^ v {do^i — cAp)\ = 0, or, explicitly, 

d^d^o, + □ fa - ci) = 0. (C30) 

Note here that under an infinitesimal U(l) gauge transformation, 5Am — * SAm + c*m£ and 8ai — * Soi + eE, where 
S is independent of <f>, the above equations of motion (and boundary conditions in the main text) are invariant. 



APPENDIX D: ON C AND F w 

The continuity of the gauge field imposes 



A 

From the (<fi(J,) component of the Israel conditions we obtain 

2fi 



0. 



¥1 



e£M 2 aff l 



)ai - cAj 



= 0, 



(Dl) 



(D2) 
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where <2j := a(w,) and fa := f(wi). From the component of the Maxwell junction conditions we get 

lev 2 



" V l fx 
1/5 { S °i 



eA 



We can eliminate (6<Ji — eA)\ iBi from Eqs. (|D2|) and (|D3|) to get 

2 ^ 



tB ( (eviiM) 2 



0. 



(D3) 



(D4) 



Equations (fC25|) . ([C26]) . (fC28|) . (|C29|> . and jD4} supplemented with appropriate boundary conditions at the poles 
form a closed set of equations that can determine C and F w . Then Eqs. (|C30p . (|D1 1) . and (|D2|) [or (|D3[) ] arc enough 
to determine Jcr^ and ^4 on the brane. In Appendix [D] we show that C and F w modes vanish everywhere both for the 
single brane and two-brane models. 

Equation (|C25[) shows that only the zero mode is present for C. The same is true for F w , as is clear from Eq. (|C29[) . 
The general solution to Eqs. (|C26|) and (|C28|) is given by 



C = 
F = 

■*- m 



C f (z) 

IM 2 F 1 (a:) 
2fi a 2 f ' 



F 7 (a;) / 



(D5) 
(D6) 



where C 1 and F 1 are integration constants. The regularity at the poles requires that and F ± vanish. Thus in the 
single brane model C = and F w = everywhere. In the two-brane model, we can write 



C = ^ +F o (a .) [ W+ J-dw, 



a 2 P 
tM 2 F°(x) 



2fi a 2 } ' 

for u>_ < w < to+. However, the junction condition (jD4| at both branes reads 



C° + 



dw I a 2 / 2 + 



^(ew+) 2 a 4 

A)-/?- / 1/2 
^(eu_) 2 a 4 



F° = 0, 



0, 



(D7) 
(D8) 

(D9) 
(D10) 



leading to C° = and F° = 0. Thus we see that (7 = and F w = everywhere also in the two-brane model. Then 
from Eq. (|E>2[1 or Eq. (|D3|) we have Sai — eA = 0, at the branes. Eq. (|C30[) now reduces to d^d^Sai = 0, which is 



solved by boi — Sa-\ 1 \x) + daf 1 {x)4>. However, Sa^' — because of periodicity in <p of Ui and the fact that 8<Ji is 
infinitesimal. Thus we get 5(Ji = eA^i = 5a^{x). 
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